OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Induction of the CLOCK Gene by E2-ERa Signaling fjfi 
Promotes the Proliferation of Breast Cancer Cells cros^rk 

Liyun Xiao^ ^ Alan K. Chang^ ^ Ming-Xi Zang^^ Hailian Bi\ Shujing Li\ Miao Wang\ Xinrong Xing\ 
Huijian Wu^'^* 

1 School of Life Science and Biotechnology, Dalian University of Technology, Dalian, China, 2 School of Life Science and Medicine, Dalian University of Technology, Panjin, 
China, 3 Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Zhengzhou University, Zhengzhou, China 



Abstract 

Growing genetic and epidemiological evidence suggests a direct connection between the disruption of circadian rliythnn 
and breast cancer. IVloreover, the expression of several molecular components constituting the circadian clock machinery 
has been found to be modulated by estrogen-estrogen receptor a (E2-ERa) signaling in ERa-positive breast cancer cells. In 
this study, we investigated the regulation of CLOCK expression by ERa and its roles in cell proliferation. 
Immunohistochemical analysis of human breast tumor samples revealed high expression of CLOCK in ERa-positive breast 
tumor samples. Subsequent experiments using ERa-positive human breast cancer cell lines showed that both protein and 
mRNA levels of CLOCK were up-regulated by E2 and ERa. In these cells, E2 promoted the binding of ERa to the EREs 
(estrogen-response elements) of CLOCK promoter, thereby up-regulating the transcription of CLOCK. Knockdown of CLOCK 
attenuated cell proliferation in ERa-positive breast cancer cells. Taken together, these results demonstrated that CLOCK 
could be an important gene that mediates cell proliferation in breast cancer cells. 
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Introduction 

Breast cancer is one of the most prevalent causes of cancer 
death among women. Prolonged exposure to estrogen is thought 
to be a major factor contributing to the development and 
progression of breast cancer [1,2]. About 70% of breast cancers 
are estrogen-dependent. Moreover, clinical studies in which anti- 
estrogen or aromatase inhibitors are used to decrease the rate of 
local and distant relapse have demonstrated that estrogen can 
facilitate the progression of breast cancer [3] . 

The molecular mechanism of breast cancer induced by estrogen 
is thought to occur through the binding of estrogen to the 
transcription factor estrogen receptors (ERs), which then binds to 
estrogen response elements (EREs) in the promoters or regulatory 
regions of target genes. ERs contain two isoforms, ERa and ERP, 
and each is encoded by a different gene. ERa is highly expressed 
in ER-positive breast cancer and is associated with breast cancer 
growth [1,4,5]. ERP is also expressed in breast cancer, but its role 
is still elusive [6]. Moreover, ERa can bind to the promoter or 
regulatory regions of target genes that contain imperfect or 
truncated EREs, and activate their transcription [7,8]. E2-ERa 
signaling plays a critical role in cell proliferation [9] . E2 promotes 
the proliferation of breast cancer cells through a number of 
estabhshed pathways [3] . 

Circadian rhythm is conserved across a wide range of 
organisms, including Arabidopsis, Drosophila, and mammals [10]. 



The duration of a circadian cycle is about 24 h. In mammals most 
physiological and behavioral functions are influenced by circadian 
rhythm. These rhythms are directed by endogenous clocks residing 
in the hypothalamic suprachiasmatic nucleus (SCN) and periph- 
eral tissues [11,12]. The molecular components of the circadian 
rhythm, clock genes and their products form the transcription- 
translation feedback loops. Two core transcription factors, 
CLOCK (circadian locomotors output cycles kaput) and BMALl 
(brain and muscle ARNT (aryl hydrocarbon receptor nuclear 
translocator)-like protein 1), form a heterodimer that binds to the 
E-box in the promoters of their target genes and activate the 
expression of these genes, including Period [PERI, 2 and 3) and 
Cryptochrome [CRTl and 2). PER and CRY proteins can form 
heterodimer complexes that translocate to the nucleus, where they 
interferes with the transcriptional activity of BMAL 1 / CLOCK to 
limit their own expression, thereby constituting a negative- 
feedback loop [13,14]. The circadian negative-feedback loop 
results in the circadian expression of clock genes. 

There is accumulating evidence to suggest that circulating 
hormones could regulate the circadian oscillations of clock gene 
expression in some brain regions and peripheral tissues [15-18]. 
The ability of rhythmically-produced hormones to regulate the 
expression of clock genes in specific tissues implies a relationship 
between circadian clock and hormone production [12,18]. 
Circulating levels of hormones can modulate circadian clock, 
which in turn regulates the periodic release of these hormones 
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[19]. Recent reports have suggested that the circadian rhythm and 
the physiological condition of the body can mutually influence 
each other in mammals [20]. More and more evidence is 

suggesting that circadian disruption is associated with tumor 
occurrent:e, including breast cancer [21-23]. Estrogen plays a 
critical role in normal mammary gland physiology. At the same 
time it also acts as a potent mammary mitogen. Although the 
circadian clock is linked to the activity of estrogen, the molecular 
mechanisms underlying the regulation of the core clock genes that 
regulate the mammary circadian regulation are largely unknown. 

It has been reported that upon treatment with E2, expression of 
the core clock gene BMALl expression is enhanced [16]. Another 
circadian clock gene, PER2, is also a target gene of ERa that is 
regulated by E2 [18,24,25]. In rat uterus, E2 induces the high 
expression of Perl [15,26]. These studies indicate that there is a 
relationship between E2-ERa signaUng and the gears of circadian 
rhythm machinery. In our previous study, we have confirmed that 
CLOCK interacts with ERa and enhances its transcriptional 
activity [27]. Another breast cancer-associated protein, DECl, has 
been shown to repress the transcriptional activity of CLOCK [28]. 
Many studies have focused on the transcriptional activity of 
CLOCK, but the transcriptional regulation of CLOCK is largely 
unknown until the revelation that the nuclear receptor REV- 
ERBa, a critical component of the circadian clock [29], is a 
transcriptional repressor of CLOCK [29,30]. Based on the 
relationship between E2-ERa signaling and circadian clock genes, 
we wanted to know whether CLOCK is under the control of E2- 
ERa signaling. 

In this study, we found that in ERa-positive breast cancer cells, 
E2 treatment increased while knockdown of ERa decreased the 
expression of CLOCK. In addition, we showed that ERa could bind 
to CLOCKvia EREs and activate CLOCA' transcription in response 
to E2. Taken together, our data suggested that CLOCK is a 
transcriptional target of ERa, and that the product of this gene can 
modulate cell proliferation in ERa-positive breast cancer cells. 

Material and Methods 

Ethics Statement 

AU research involving human participants have been approved 
by the institutional review board of Qiqihar Medical University. 
We have obtained the written consent from all the human 

participants and our clinical investigations have been conducted 
according to the principles expressed in the Declaration of 
Helsinki. 

Cell Culture and Experiment Reagents 

Human breast cancer cell lines MCF-7, T47D, MCFIOA and 
MDA-MB-231 have been used in our previous studies [31-34]. 
MCF-7 and MDA-MB-231 cells were cultured in Dulbecco's 
Modified Eagle's Medium (DMEM) containing 10% fetal bovine 
serum (Hyclone, Logan, UT), 100 Hg/ml penicillin and 100 |Xg/ 
mi streptomycin. T47D ('("lis were cultured as pre\'iously described 
[35]. MCFIOA were cultured in DMEM/F12 (1:1) containing 5% 
horse serum (Hyclone, Logan, UT), 100 ng/ml Cholera Toxin, 
10 ng/ml bovine insulin, 0.5 |J,g/ml hydrocortisones (sigma) and 
20 ng/ml EOF (Peprotech, Rehovot Israel). Unless otherwise 
stated, all cell cultures were incubated at 37°C in the presence of 
5% CO2. Cycloheximide (CHX), actinomycin D (Act D) and anti- 
FLAG antibody were obtained from Sigma. Anti-ERa, anti-FIA, 
and anti-CLOCK were obtained from Santa Cruz Biotechnology 
(Santa Cruz, CA, USA), and 17(3-estrogen (E2) and ICI 182780 
(ICI) were obtained from Abeam (Cambridge, UK). 



Plasmlds Constructions 

The promoter region of human CLOCK (gene ID 9575) was 
amplified from a human genomic DNA by PGR and cloned into 
the plasmid pGL3-basic. Two truncated versions of CLOCK 
promoter were constructed, and each was fused to a luciferase 
reporter gene. CLOCK-WT-Luc (-884/-I-992) was amplified by 
sense primer 5'-GATCGGTACCCCAGTAGAAGCACT- 
GAAATG-3' and antisense primer 5'-GATCCTC- 
GAGTCGCTGGAGTCAGACGCTAAT-3'; truncated 
CLOCK-Ml-Luc (— 297/+63) was amplified by sense primer 5'- 
GATCGGTACCAAAGCCAAAGAGCCTCC-3' and antisense 
primer 5'-GATCCTCGAGTTTTAAACCGGCAGCC-3'; trun- 
cated CLOCK-M2-Luc (-I-552/-I-992) was amplified by sense 
primer 5'-GTACGGTACCGAGCTGCGGCCGATTCC-3' and 
antisense primer 5'-GATCCTGGAGTCGCTGGAGTCA- 
GACGCTAAT-3'. An ERE half-site and an ERE in CLOCK- 
WT-Luc were also mutated using a site-directed mutagenesis kit 
(Stratagene, La JoUa, CA, USA) according to the manufacturer's 
instruction. The mutant CLOCK-M3-Luc contained base substi- 
tutions in the ERE half-site of CLOCK (-I-2/+6), whereas the 
mutant CLOCK-M4-Luc contained base substitutions in the ERE 
of CLOCK (+753/+764), while the mutant CLOCK-M5-Luc 
contained base substitutions in both ERE half-site and ERE. 
The mutants were each generated by a pair of primer: CLOCK- 
M3-Luc (sense primer 5'-CCGCGGGGTCGCTTGCGACG- 
CATGCGCCGG-3', and antisense primer 5'- 
CCGGCGCATGCGTCGCAAGCGACCCCGCGG-3'), 
CLOCK-M4-LUC (sense primer 5'-CTGGGGAGCCGCTAGG- 
CAATGTTGCGCACTTTATTCCTGTCA-3', antisense prim- 
er 5'-TGACAGGAATAAAGTGCGCAACATTGCC- 
TAGCGGGTCCCCAG-3'). AU cloned and mutated genes were 
verified by DNA sequencing. HA-REV-ERBa construct was a gift 
kindly provided by Dr. Hiromitsu Negoro (Kyoto University 
Graduate School of Medicine). 

Luciferase Reporter Assay 

HeLa or MCF-7 cells were transfected with the appropriate 
plasmids, and 24 h after transfection, the cells were rinsed with 
PBS and subjected to luciferase activity assays. Briefly, the cells 
were lysed in cold buffer containing 25 mM glycylglycine 
(pH 7.8), 1% Triton X-100, 4 mM EGTA, 1 mM DTT, and 
15 mM MgS04. Five microliters of assay bufier (1 M MgClj, 
0.5 M KH2PO4, and 0.1 M ATP) and 100 |.il 0.2 mM luciferin 
potassium salt (BD Biosciences Pharmingen, Franklin Lakes, NJ) 
were added to 45 ^ll cell lysate [36], and the luciferase activity of 
the sample was measured with a Centro LB 960 Microplate 
Luminometer (Berthold Techologies GmbH Co KG, Germany). 
The efficiency of transfection was evaluated by transfecting the 
cells with a fi-galactosidase construct. Briefly, 20 |J.l cell lysate was 
added to 50 ^1 (5-galactosidase bufffer (60 mM Na2HP04-12H20, 
40 mM NaH2P04-2H20, 10 mM KCl, 1 mM MgS04-7H20, 
and 6 mg/ml ONPG) and the absorbance of the sample was 
measured at 450 nm. 

RNA Extract and RT-PCR 

Triozol reagent (Invitrogen, Grand Island, NY) was used to 
extract the total RNA from MCF-7 cells. The extraction was 
performed according to the manufacturer's instruction and the 
concentration of RNA was quantified by optical density. One 
microgram of total RNA was retrotranscribed into cDNA using 
Reverse Transcription System (TAKARA, Dalian, China). Real- 
time PGR was performed with a Roter-Gene 3000 (Corbett 
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Research, Australia) using the following primers as reported: 
5'-AAGTTAGGGCTGAAAGACGACGA-3' (sense) and 5'- 
GAACTCCGAGAAGAGGCAGAAG-3' (antisense) for CLOCK 
[16]; 5'-GAAGGTGAAGGTCGGAGTC-3' (sense) and 5'- 
GAAGATGGTGATGGGATTTC-3' (antisense) for GAPDH 
[16]. The cDNA was eombined with the appropriate pair of 
primers and the maxima SYBR Green qPCR Master Mix 
(Thermo scientific) and subjected to the following reaction: initial 
denaturation step of 95°C for 10 min; and 40 cycles of 95°C for 20 
s; 56°C for 20 s and 72°C for 20 s. The eEBciency of the real-time 
PGR assay was determined from the amplification efficiency E and 
linear correlation coefficient R^. Ten-fold serial dilutions (from 
10"'' to 10"^) of cDNA generated from MCF-7 cells were used in 
the real-time PGR assay to generate a set of data for the standard 
curve. E and R^ values were calculated from the standard curve as 
in previous report [37]. To evaluate the quality of the product of 
real-time PGR, melt curve analyses were performed after each 
reaction. GAPDH is a frequently used housekeeping gene in real- 
time PGR as it is expressed at a relatively constant level in various 
tissues, including breast tissue, under normal and pathophysiolog- 
ical conditions [38-40]. Therefore, the expression level of CLOCK 
was normalized the expression level of GAPDH using Roter Gene 
6.0 software. Relative expression was determined using the 2 
method with GAPDH as the reference gene. Each target was 
measured in triplicate. 

Western Blot Analysis 

Preparation of cell extracts and subsequent western blot analysis 
were carried out as previously described [28]. Immunoblot data 
were quantified by scanning the appropriate bands of interest and 
plotted as relative density of gray scale. 

RNA Interference 

ERa shRNA-expression vector was constructed by DNA vector- 
based shRNA synthesis using the vector pRNATU6. 1 (GenScript, 
Piscataway, NJ). The sequence of ERa used for knockdown study 
were 5'-GCTACTGTTTGCTCCTAAC-3' (shERQ(#l) [41] and 
5'-AGTTTGTGTGCCTCAAATC-3' (shERa#2) [42]. The 
sequences used for silencing the expression of CLOCK have been 
described in our previous study [27]; and the sequenc:e of the 
control shRNA is 5'-GACGCTTACCGATTCAGAA-3' [35], 
which has no significant homology with human gene sequence. 
shERa# 1 and shERa#2 expression vectors were verified by DNA 
sequencing. 

Chromatin Immunoprecipitation Assays 

MCF-7 cells were grown for 2 days in phenol red-free DMEM 
containing 5% charcoal-dextran-stripped FBS. The cells were then 
treated with or without 1 |iM E2 for 1 h, and then crosslinked with 
1% formaldehyde in PBS for 1,5 min at room temperature. Crude 
cell lysate was sonicated to generate DNA fragments of 300 to 
1500 bp. The generated DNA fragments were diluted 1:10 in 
dilution buffer (150 mM NaCl, 2 mM EDTA, 1% Triton X-100, 
and 20 mM Tris-HCl pH 8.0) [35]. Protein A and anti-ERa 
antibody or rabbit IgG were then added to the diluted sheared 
chromatin, and the mixture was incubated with constant rotation 
at 4°C for overnight. The immunoprecipitated chromatin was 
purified from the chromatin-antibody mixture and eluted in 
elution buffer (50 mM Tris-HCl pH 8.0, 10 mM EDTA, and 10% 
SDS). The isolated DNA was subjected to PGR to amplify the 
regions using specific primers: 5'-GAGCTGCGGCCGATTCC- 
3' (sense) and 5'-GCTGCTCCAAACGTGC-3' (antisense) for 
CLOCK (+672/+805); 5'-AAAGCCAAAGAGCCTCC-3' (sense) 
and 5'-TTTTAAAGCGGCAGCC-3' (antisense) for CLOCK 



(-297/-I-63); and 5'-TGAAAGAGGGAGGAGTCAAAGAT-3' 
(sense) and 5'-AGCAAGACGGAGGCAAAGTTATT-3' (anti- 
sense) for CZOCr(-1866/-1626). Total input DNA (1:10 dilution) 
was used as a positive control for the PGR reaction. The anti-IgG 
antibody was used as a non-specific binding control. The PCR 
products were analyzed by electrophoresis using 1.5% agarose gel. 

Cell Proliferation Assays 

Cell proliferation was assessed by MTT assay. MCF-7 or T47D 
cells were cultured for 24 h in phenol red-free DMEM 
supplemented with 5% charcoal-dextran-treated fetal bovine 
serum. The cells were then transfected with shControl (shCon), 
shCLOCK or shERa#l construct. After 24 h, the cells were 
treated with vehicle or 1 ^M E2 for several days, and then 
subjected to MTT assay performed with a commercial kit (Key 
Gen) according to the manufacturer's protocol. TIk; absorbance of 
the samples was read at 490 nm [27]. For colony formation assays, 
MCF-7 cells were transfected witii shControl (shCon), shCLOCK 
or shERa#l construct, and the cells were then collected and 
plated at a density of 1000 cells/well in 24-well plates, and treated 
with \'ehicle (ethanol) or 1 |iM E2 for seven days. After that the 
cells were washed with PBS, fixed with ethanol, stained with 0. 1 % 
crysteJ violet, and then photographed. The stained cells were 
solubilized in 10% SDS, and absorbance was measured at 570 nm 
[43]. 

Soft-Agar Colony Culture 

The anchorage-independent growth of MCF-7 cells was 
estimated by soft-agar colony culture as described previously 
[44,45]. MCF-7 cells were transfected with CLOCK expression 
vector, ERa cx])rcssion vector or empty vector (pcDNA3) and 
grown in the presence of 1 |J,g/ml G418 for 2 weeks. Aliquot of 
the cell suspension containing 1000 cells was mixed with 1 ml 
DMEM containing 10% fetal bovine serum and 0.4% agar, and 
then poured over a layer of solidified 0.7% agar (prepared in 1 ml 
medium) in a well of a 6-well plate. Additional 500 |J.l DMEM 
containing 10% fetal bovine serum was added to the well every 
two days. One week after seeding, photographs of the colonies 
were taken under phase-contrast microcopy, and the diameters of 
the colonies were measured by the software Image-Pro Plus 6.0 
[46]. 

Breast tumor tissue samples for immunohistochemical 
assay 

Breast tumor samples used for immunohistochemical assay were 
obtained from Qiqiliar Mc'dical University. The spc'cimens were 
obtained from female patients of Han Chinese descent, with ages 
ranging from 39 to 75 years old (average age of 56.7 years). A total 
of 32 specimens were obtained and 19 of these were ERa-positive 
and 13 were ERa-negative, as determined by clinical diagnosis 
performed by Qjqihar Medical University. The tumor grades were 
recorded as II (22 specimens). III (3 specimens) or II-III for 
obscure tumor grade (7 specimens). 

Immunohistochemical Assay 

All the obtained human breast tumor specimens were analyzed 
by immunohistochemical assay, which was performed with the aid 
of an immunohistochemical assay kit (Maixin Bio, China). Sections 
of the tissues were first fixed in 10% buffered formalin. The fixed 
tissues were embedded in paraffin and then deparafinninized 
before being rehydrated using standard procedures [47]. The 
endogenous peroxidases of the samples were quenched with H2O2 
in methanol. After that, the samples were incubated in blocking 
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solution (4% bovine serum albumin, 0.1% Triton X-100, 0.1 M 
PBS) for 10 min [48] and then stained with citrate buffer (pH 6.0) 
containing a 1:100 dilution of either rabbit anti-human CLOCK 
(Santa Cruz, CA, USA) or rabbit anti-human ERa (Santa Cruz, 
CA, USA) for overnight at 4°C. Following incubation with 
primary antibody, the slides were then incubated with a 1:200 
dilution of biotinylated goat anti-rabbit IgG for 10 minutes at 
room temperature, and then incubated in Avidin-Biotin Peroxi- 
dase Complex for 10 minutes at room temperature. After washing 
for three times in PBS they were incubated with DAB 
(diaminobenzidine), which was used as chromagc'n for the 
antibody. The intensity of the staining, which reflected the level 
of CLOCK/ERot in the sample, was quantified by assigning it an 
H score. The procedure for acquiring an H score was performed 
as described previously [49,50] . The levels of CLOCK expression 
in these breast tumor samples were therefore sorted according to 
their H scores. 

Statistical Analysis 

A Chi-square (y^) test was used to examine the correlation 
between CLOCK and ERot gene expression in breast cancer 
tissues from 32 patients. All other data were expressed as means ± 
SDs. Differences between mean values were analysed by ANOVA, 
followed by the Bonferroni test for pairwise comparisons. 
Statistical significance was considered at the P<0.05 level. 

Results 

CLOCK protein is upregulated in ERa-positive breast 
tumor 

It has been reported that aberrant ERa signaling is related to 
the occurrence of ERot-positive breast tumor. ERa-positive breast 
tumor generally has a better prognosis, and is responsive to anti- 
estrogen therapy. However, the role of CLOCK in breast tumor 
has not been elucidated. In order to examine the relationship 
between CLOCK and ERa in breast tumor, we compared the 
protein levels of CLOCK and ERa in ERa-positive breast tumor 
samples with those of ERa-negative breast tumor samples (Fig. 1 A). 
A total of 32 tissue samples (19 ERa-positive and 13 ERa- 
negative) were analyzed by immunohistochemical assay. Fourteen 
of the ERa-positive samples showed high CLOCK expression 
(74%). As for the 13 ERa-negative samples, high CLOCK 
expression was found in only 6 samples (46%) (Fig. IB). The data 
appeared to suggest a correlation between the ERa and CLOCK 
in ERa-positive breast tumors. 

ERa regulates the level of CLOCK protein 

Circulating estrogen (E2) is known to modulate the expression 
of some clock genes. Given the apparent correlation between ERa 
and CLOCK expression in the breast-tumor samples analyzed, we 
next examined whether CLOCK expression could be stimulated 
by E2. Treatment of the ERa-positive breast cancer cell lines, 
MCF-7 and T47D ceUs with 1 nM E2 for 24 h resulted in 
increased expression of CLOCK protein in these cells, but the 
same treatment given to the ERa-negative cell lines, MDA-MB- 
231 and MCFIOA resulted in no apparent effect on CLOCK 
expression (Figs. 2Ac&B). The importance of ERa with respect to 
increased expression of CLOCK when MCF-7 and T47D cells 
were treated with E2 was further investigated by observing the 
change in CLOCK expression when these cells were treated with 
tiie anti-estrogen agent, ICI182780 (ICI) instead of E2. ICI 
treatment resulted in a reduction of CLOCK expression. 
Furthermore, the E2-enhanced expression of CLOCK in these 
cells was partiy reversed when the same cells were also treated with 



ICI (Fig. 2C). The level of ERa in these cells was also down- 
regulated after treatment with E2 or ICI, which was consistent 
with previous reports [51]. In contrast, MDA-MB-231 cells, which 
are ERa-negative but ERpi-positive, showed no obvious changes in 
the level of CLOCK expression when treated with E2 or ICI 
(Fig. 2C). T47D is an ERa-positive but ERfi-negative breast 
cancer cell line. Overexpression of ERP in T47D cells had no 
profound impact on the level of CLOCK (Fig. 2D). These results 
suggested that ERP may have a minimal effect on the expression 
of CLOCK. To determine whether E2 could stimulate the 
expression of CLOCK through ERa, the expression of ERa was 
knocked down with two different shRNAs, shERa^l and 
shERa#2, which target different regions of the ERa mRNA to 
avoid possible off-target effect. The effectiveness of the two 
shRNAs has already been demonstrated by other investigators 
[41,42]. Knockdown of ERa decreased the expression of CLOCK 
(Fig. 2E). When ERa was overexpressed in MCF-7 cells, the level 
of CLOCK protein increased slightly (Fig. 2F). Taken together, 
the results showed that expression of CLOCK is subject to control 
by ERa. 

ERa modulates the transcription of CLOCK 

The effect of ERa on the level of CLOCK expression was 
further investigated by determining the changes in the level of 
CLOCKmRNA in MCF-7 cells in response to E2 or ICI after 24 h 
of treatment. CLOCKmRNA level was up-regulated in response to 
E2 in a dose-dependent manner in the range of 10"'° to IO'^'M 
(Fig. 3A). Thus 10~^ M E2 was chosen for subserjuent studies. 
CLOCKmRKA level was increased 4 h after E2 treatment (Fig. 3B). 
However, when the cells were treated with ICI, the level of CLOCK 
mRNA was reduced compared to that of the control (Fig. 3C). The 
effect of ERa on the modulation of CLOCK transcription in 
response to E2 was further confirmed by overexpressing ERa in 
MCF-7 cells and knocking down ERa with shERa. As expected, 
ERa c-ctopic expression up-regulated CLOCA" transcription, while 
knockdown of ERa down-regulated CLOCff transcription (Fig. 3D). 
In addition, Act D repressed the basal expression of CLOCK and 
abolished E2-induced up-regulation of CLOCK. Although Act D 
globally represses gene transcription, including the transcription of 
GAPDII, it inhibited the transcription of CLOCK more, than that of 
GAPDH. In contrast, cycloheximide (CHX, a translation inhibitor) 
had no effect on E2-induced up-regulation of CLOCK (Fig. 3E), 
indicating that CLOCK is a primary ERa transcriptional target 
because the effect of E2 does not require the synthesis of new 
proteins since all necessary factors are preexisting in the cells. 

ERoc regulates CLOCK promoter activity 

The transcription of CLOCK has been shown to be repressed 
by nuclear receptor REV-ERBa through its interaction with the 
REV-ERB response element (RevRE) located 760-771 bases 
downstream the transcription start site (TSS) of the CLOCK gene 
[30]. As the location of the regulatory element is in the first intron 
of the CLOCKgene, a 1877-bp fragment encompassing the regions 
upstream and downstream of the TSS was cloned (See Fig. 4A). 
This fragment which included the reported RevRE was fiised to 
the luciferase g(;nc. gciu-rating the construct CLOCK-WT-Luc. 
HeLa cells transfixtcd with CLOCK-WT-Luc and REV-ERBa 
showed a decreasing trend in luciferase activity that was 
dependent on the dosage of REV-ERBa (Fig. 4B), consistent with 
former report. However, overexpression of ERa could attenuate 
the repression of luciferase activity by REV-ERBa (Fig. 4C). HeLa 
cells were used because these [:ells have no detectable levels of ERa 
or ERp. The cells were transfected with increasing amounts of 
ERa-expression plasmids and CLOCK-WT-Luc exhibited en- 
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Figure 1. Correlation between ERa and CLOCK expression in human breast tumor tissue samples. A, Representative results showing the 
immunohistochemical staining of ERa and CLOCK in serial sections of the breast tumor tissues. Each sample was incubated with antibody against ERot 
or CLOCK. Positive staining and negative staining are indicated by brown and blue staining, respectively (x200 Magnification). B, Correlation between 
ERa and CLOCK expression suggested by the 32 breast tumor samples, test was used for statistical analysis. P values less than 0.05 were considered 
to indicate statistical significance. 
doi:1 0.1 371 /journal.pone.0095878.g001 



hanced expression of luciferase activity that paralleled with the 
dosage of ERot gene (Fig. 4D). HeLa cells transfected with 
CLOCK-WT-Luc plus increasing amounts of ERP resulted in 
minimal increase in luciferase activity compared to cells transfect- 
ed with CLOCK-WT-Luc and ERa (Fig. 4E). To determine 
whether ERa could stimulate CLOCA" transcription, the luciferase 
activity of MCF-7 cells transfected with CLOCK-WT-Luc 
followed by treatment with E2 or ICI was determined. The level 
of CLOCK-WT-Luc luciferase activity was increased by 3.8-fold 
in the presence of E2, but such an enhancement was attenuated by 
ICI (Fig. 4F). In another breast cancer cell line T47D, similar 
results were obtained (Fig. 4G), but in MDA-MB-231 cells, the 
level of luciferase activity expressed by CLOCK-WT-Luc did not 
change significantly in response to E2 or ICI (Fig. 4H). As for 
MCF-7 cells, knockdown of ERa with shERa resulted in 
decreased level of reporter activity (Fig. 41). These results 
confirmed that CLOCK was transcriptionally regulated by ERa. 



ERa bounds to CLOCK promoter regions in response to 
E2 

To map the ERa responsive regions within the CLOCK 
promoter region, a computational analysis was performed and 
the results indicated that half estrogen response element (V2ERE) 
was present at the +2 to +6 region while an ERE was present at the 
+753 to +764 region. To delineate which portion of the promoter 
was responsive to ERa, two truncated versions of the promoter- 
fused luciferase were constructed, CLOCK-Ml-Luc (-297/+63) 
and CLOCK-M2-Luc (+552/+992), and their activity in response 
to ERa was tested. MCF-7 cells transfected with either construct 
showed a decreased level of luciferase activity when ERa was 
knocked down, a trend that was also exhibited by wide-type 
CLOCK-WT-Luc (Fig. 5A). In a different experiment, constructs 
of CLOCK luciferase reporter bearing point mutation in the 
CLOCK component were also made. The point mutation 
consisted of two nucleotide substitutions at either the half ERE 
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Figure 2. Western blot analyses of CLOCK and ERo expression in cells treated with E2 or ICI. A, ERoc and ERp expression in MCF-?, T47D, 
MDA-MB-23 1, and MCF10A cells. B, CLOCK expression in MCF-7,T47D, MDA-MB-23 1, and MCF10A cells that had been treated with vehicle (control) or 
1 |iM E2 for 24 h. Cells were cultured in 5% charcoal striped FCS and phenol red free medium for 2 days before stimulated with E2. C, CLOCK and ERot 
expression in MCF-7, T47D and MDA-MB-231 cells that had been treated with vehicle, 1 |iM E2 or 0.1 |j.M ICI alone or in combination for 24 h. Cells 
were cultured for 2 days in 5% charcoal striped FCS and phenol red free medium for two days before they were treated with ER ligands. D, CLOCK and 
ERa expression in T47D cells transfected with empty vector pcDNA3 or pcDNA3-Flag-ERp. E, CLOCK and ERa expression in MCF-7 cells transfected 
with control shCon or two different shERa (shERa#1 and shERa#2). F, CLOCK expression in MCF-7 cells transfected with pcDNAS or pcDNA3-Flag- 
ERa. B-F, 24 h after transfection, the cells were harvested and subjected to western blot analysis. In all experiments (A-F), (J-actin expression was used 
as a reference. The blot shown is the representative result from three independent experiments. Image of the blot is shown in the top panel of each 
figure, with the quantitative analysis of the bands in the blot shown in the plot below. The levels of CLOCK or ERa signal obtained from control cells 
were set to 1 . All experiments were repeated at least three times. Data shown in the graphs are the means ± SDs of three experiments. P value was 
determined by ANOVA with Bonferroni test (*, P<0.05. ns, not significant). 
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(CLOCK-M3-Luc), six nucleotide substitutions at the ERE site 
(CLOCK-M4-Luc) or both (CLOCK-M5-Luc). HeLa ceUs 
transfected with the wild-type construct and those transfected 

with any of the three mutant forms showed similar levels of 
luciferase activity in the absence of ERa overexpression. With 
ERa overexpression, the level of luciferase activity increased by 
about eight fold in the case of wild-type construct, about six fold 
for CLOCK-M3-Luc, four fold for CLOCK-M4-Luc and three 
fold for CLOCK-M5-Luc (Fig. 5B), indicating that although the 
presence of intact ERE was important for CLOCK-driven 
luciferase activity, such activity was dependent on the presence 
of ERa, which interacted with the ERE of the CLOCK promoter 
region. It was worth noting that ERa still activated the activity of 
CLOCK-M5-Luc despite the absence of any ERE (Fig. 5B). This 
could possibly be due to ERa interacting with other transcription 
factors that could bind to other regions of the CZOCff promoter in 



CLOCK-M5-Luc [52,53]. Such binding would effectively enable 
ERa to bind to CLOCK-M5-Luc indirectly and consequentiy 
activate the transcription of the reporter gene. Surprisingly, in the 
case of nucleotide-substitution mutation, the result appeared to 
indicate that interaction between ERa and ERE was more 
important than between ERa and the half ERE, whereas in the 
case of truncation mutation, the result seemed to indicate the 
opposite. Nevertheless, both did confirm that ERE was essential 
for ERa-mediated upregulation of CLOCK activity. 

We further tested the binding of ERa on the CLOCK gene in 
MCF-7 cells by chromatin immunoprecipitation (ChIP) assay in 
vivo. Following treatment of MCF-7 cells with E2, the DNA of the 
cells was immunoprecipitated and analyzed by PCR using the 
CLOCK gene, specific primers correspond to regions —297 to +63, 
+672 to -H805, or -1866 to -1626. ERa only bound to regions 
-297 to -F63 and -F672 to -1-805, but not region - 1866 to - 1626 
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Figure 3. ERa ligands regulate the expression of CLOCK at the transcription level. Analysis of CLOCK mRNA levels in l\/\CF-7 cells by real- 
time PCR. IMCF-? cells were cultured in phenol red free medium and charcoal striped FCS medium for 2 days before being treated with E2 or ICI and 
the expression of CLOC/f was then analyzed by real-time PCR. Expression of CLOC/f was normalized against GAPDH mRNA level (internal control). A, 
Cells treated with different concentrations of E2 (10"'° to 10"® IVI) for 8 h. B, Cells treated with 1 nM E2 for different periods of time. C, Cells treated 
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repeated at least three times. Data shown are the means ± SDs. P value was determined by ANOVA with Bonferroni test (*, P<0.05. ns, not 
significant). 
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in the presence of E2 (Fig. 5C). Regions —297 to H-63 and -H672 to 
-1-805 contains the VjERE and ERE, respectively (left panel of 
Fig. 5A). These results suggested that the binding of endogenous 



ERot to the CLOCK gtne is dependent on E2. ChIP results further 
supported the involvement of ERa in the transcription of CLOCK 
gene. 
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Figure 4. ERa regulates CLOCK promoter activity. A, Schematic illustration of estrogen response elements in CLOCK promoter containing CLOCK 
sequence from -884 to 4-992 fused to luciferase (CLOCK-WT-Luc). B, Luciferase activity of HeLa cells transfected with CLOCK-WT-Luc and increasing 
amounts of REV-ERBa expression plasmid. C, Luciferase activity of HeLa cells transfected with CLOCK-WT-Luc plus REV-ERBa or ERa expression 
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CLOCK is required for the proliferation of breast cancer 
cells 

The potential relevance of our lindings to the biology of breast 
cancer cells was investigated by looking at the effect of reduced 
CLOCK expression on cell proliferation. E2 treatment stimulated 
the proliferation of MCF-7 cells transfected with shControl 
(Fig. 6A). In contrast, knockdown of CLOCK or ERa inhibited 
cell proliferation in the absence and presence of E2 (Fig. 6A). 
Similar results were obtained for T47D cells (Fig. 6B). These 
results suggested that CLOCK and ERot could promote cell 
proliferation. We next examined the effect of CLOCK knockdown 
on the colony formation of MCF-7 cells. E2 treatment increased 
the colony formation of MCF-7 cells (top panel of Fig. 6C). 
Knockdown of CLOCK decreased the colony formation of MCF- 
7 cells in the absence and presence of E2 compared to shControl 



(Fig. 6C). Knockdown of ERa expression also decreased the 
colony formation ability of the cells as expected, and at the same 
time the cells lost the response to E2 (bottom panel of Fig. 6C). 
Three-dimensional cell culture is considered superior to mono- 
layer cultures because the growth of a cell in a three-dimensional 
culture resembles more the growth in an in vivo environment. Soft- 
agar culture was chosen to test the difference in growth among the 
MCF-7 cells transfected with different vectors (empty vector or 
vector harboring CLOCK or ERa gene). Cells that overexpressed 
CLOCK or ERa showed more growth than the cells that were 
transfected with the empty vector (Fig. 6D). These results indicated 
that the induction of CLOCK expression via ERa appears to 
constitute a driving force in the proliferation of ERa-positive 
breast cancer cells. 



PLOS ONE I www.plosone.org 



8 



May 2014 I Volume 9 | Issue 5 | e95878 



E2-ER!x Positively Regulates CLOCK 



A /^'AEBE /S/ ERE 



Luc 



&■ 



pGL3 

CLOCK-WT-Luc 
CLOCK-MI -Luc (-297/+63) 

— j Luc I CL0CK-M2-LUC (+552/+992) 

\'AEBE /I/ ERE j^Mnt'AERE AfMutERE 
A / gL— I Luc 




5 10 15 20 25 

Luciferase Activity (xIO^) 



4—^ 



Luc 



4 — 



Luc 



CLOCK-WT-luc 
CLOCK-M3-LUC 
CL0CK-M4-LUC 

CL0CK.-M5-LUC 



aERa 




3 6 9 

Relative Luciferase Activity 




CLOCK 



-297 



CLOCK -297 to +63 



+63 



CLOCK +672 to +805 



CLOCK -1866 to -1626 



n 



+672 +805 



-1866 



-1626 



•M: 



n 



Figure 5. ERa bounds to CiOOf promoter regions in response to E2. A, Schematic representation of the ERE sites within the QOCK promoter 
regions in the CLOCK-WT-Luc constructs. Constructs containing wild-type promoter and mutant promoters (truncation) are shown. Luciferase activity 
of IVlCF-7 cells transfected with the indicated constructs together with or without shERa#1 are shown on the right. B, CLOCK luciferase reporter 
constructs containing wild-type and mutant CLOCK promoters with point mutation in the EREs are shown, together with the luciferase activity of 
HeLa cells transfected with one of these constructs together with or without ER7.. A and B, all experiments were performed in triplicate and repeated 
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Discussion 

Apparently, women who have been exposed to artificial light at 
night for long term, or have been working with jobs that can alter 
their circadian rhythm are predisposed to brea.st cancer [54-56]. 
Growing evidence suggests that breast tumorigencsis is associated 
with the disruption of circadian clocks [23,56,57]. Thus it is 
important to investigate how estrogen signahng, which is vital to a 
number of cellular processes and the onset of breast cancer, is 
integrated with the circadian clock. CLOCK is a core transcrip- 
tion factor in the transcription-translation feedback loops of the 
machinery that regulates circadian rhythm. However, the 
molecular details of the transcriptional regulation of CLOCK 
remain largely unknown. In this study we demonstrated that 
CLOCK, similar to other circadian clock genes, is subject to 
modulation by estrogen in breast cancer cells. 

A higher percentage of ERot-positive breast tumor samples that 
we analyzed revealed a high level of CLOCK protein compared to 



ERa-negative breast tumor samples (74% versus 46%, Fig. IB), 
which suggested that the transcription of CLOCK in ERa-positive 
tumor may be upregulated. However, there appeared to be no 
correlation between the expression of CLOCK and ages. Due to 
the obscure tumor grades recorded for some of the specimens, it 
was not possible to establish a correlation between CLOCK 
expression and tumor grades. The time of tumor resection may be 
important for the study of circadian clock proteins, but in our 
study it was not thoroughly recorded. To our knowledge, 
differences in the expression levels of breast-tissue CLOCK 
protein between day and night, or changes in the expression 
levels during menstruation have not yet been reported. Although 
the number of breast tumor samples analyzed may be low, the 
correlation between CLOCK and ERa was statistically significant 
and should not be overlooked as a chance event, as subsequent 
experiments employing two ERa-positive breast carcinoma cell 
hnes, MCF-7 and T47D, revealed that the levels of CLOCK 
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Figure 6. CLOCK promotes IVICF-7 cells proliferation. Cells were transfected with control shRNA (shCon), shCLOCK or shERa#1 in the presence 
or absence of E2 for six or seven days followed by MTT assay or crystal violet staining. MTT assay of MCF-7 (A) and T47D (B) cells. The cells were 
treated with E2 for six days. C, Crystal violet staining (IVICF-7 cells). The cells were treated with E2 for seven days. Viable colonies were stained with 
0.1% crystal violet and photographed. The dye taken up by the colonies were solubilized in 10% SDS and quantified by absorbance at 570 nm. 
Representative images are shown on the left panel of Figure 6C, and the corresponding quantitative analyses are shown on the right panel. Only 
representative data from three independent experiments are shown. D, Representative colonies of each experimental group are shown. IVICF-7 cells 
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three times. A-C, Data are the means ± SDs. P value was determined by ANOVA with Bonferroni test (*, P<0.05). 
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protein and mRNA were indeed regulated by ERa (Fig. 2B). 
DifiFerences in response to E2 between MDA-MB-231 cells (ERP- 
positive/ERa-negative) and T47D cells (ERa-positive/ERP-neg- 
ative) suggested that ERP may be less important in the modulation 
of CLOCK in response to E2, compared to ERa. This was 
demonstrated by the lack of change in CLOCK protein level in 
T47D cells overexpressing ERP (Fig. 2D). These results were 
consistent with the reporter gene experiments. ICI competes with 
E2 for binding to ERa and this leads to promotion of ERa 
degradation, disruption of its localization to the nucleus and 
subsequent dimerization [51,58,59]. In our studies, E2 decreased 
the expression level of total ERa, but stimulated the expression of 
CLOCK (Fig. 2C). This result seemed paradoxical considering the 
relationship between ERa and CLOCK. In fact, ERa tends to 
accumulate in the nucleus upon E2 stimulation, and may enhance 
its own transcriptional activity. Subsequently, ERa is ubiquitinated 
and degraded through the 26S proteasome pathway [60,61]. Thus 
a possible mechanism could be that in ERa-positive breast cancer 
cells, E2 decreases the expression level of ERa, but at the same 
time, increases its nuclear translocation, resulting in an increased 
level of ERa in the nucleus and hence enhancement in ERa 
activity [51]. 

E2-ERa signaling mediates the transcription of target genes 
through classical and non-classical pathways. In the classical 
pathway, ERa binds E2, and becomes dimerized. The dimeric 
ERa then interacts with the conserved, imperfect or truncated 
EREs in the promoter or regulatory regions of the target genes to 



activate or repress their transcription. In the non-classical 
pathway, ERa modulates gene transcription through interacting 
with other transcription factors, such as AP-1, NF-kB or Spl 
[52,53]. Two putative ERa-binding sites were identified in the 
promoter of CLOCK. The sequence (TGACG) of the site located 
downstream the TSS site {+2 to -1-6) was the same as the ERE 
located in the Metastasis Associated protein 3 {MTA3) promot- 
er[62], while the sequence (AGGCCTTGTGACCC) of the other 
site (-H753 to -F764) overlapped with the site of RevRE [30]. The 
overlapping sequence is GTGACCC. The activation of CLOCK 
transcription by ERa was inhibited by the coexpression of REV- 
ERBa (Fig. 4C). This may be due to the consequence of the 
interplay between REV-ERBa and ERa, both of which competed 
for the cis-acting elements in CLOCK. The binding of ERa to 
CZOCA" promoter was confirmed in vivo by ChIP analysis (Fig. 5C). 
Whether ERa would direcdy bind to the CLOCiT promoter in vitro 
will be a subject of further study. Importandy, knocking down the 
expression of CLOCK attenuated the proliferation of MCF-7 
cells, leading to colony formation and soft-agar colony growth 
(Figs. 6A, C(&D). In addition, Brooke H. Miller et al reported that 
Clock mutation significantly inhibits the growth and proliferation of 
fibroblast cells derived from mouse embryos [14]. These results 
seem to provide evidence for a role of CLOCK in cellular 
proliferation. 

In ERa-positive human telomerase-immortalized breast epithe- 
lial cell line, the transcriptions of key clock genes, such as PERI, 
PER2, PER3, BlVIALl, CRTl, CRr2 and Rev-Erboi, and ESRl {ERa) 
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Figure 7. Proposed model showing the crosstall< between E2-ERa signaling and circadian rhythm. 

doi:10.1371/journal.pone.0095878.g007 



were found to display circadian oscillation after entrainment, 
which was applied using a serum shock method [16,63]. In 
contrast, ERot-positive and ERot-negative breast cancer cells show 
a disrupted inner clock following entrainment [16]. Moreover, 
ESRl mRNA level in ERot-positive breast cancer cells, such as 
MCF-7 and T47D does not show circadian oscillation [16]. It is 
attractive to elucidate whether the loss of circadian oscillation of 
ERo( may actually contributes to the abnormal expression of 
CLOCK and cell proliferation in breast cancer cells. Since 
CLOCK is a core transcription factor in mammalian circadian 
clock, it is reasonable to speculate that abnormal activation of E2- 
ERot signaling could induce the overexpression of CLOCK, and 
disrupt the circadian clock in breast cancer cells. It is worthwhile 
to note that the present work examined the transcriptional 
mechanism of CLOCK in ERa-positive breast cancer cells, and 
tried to determine if a correlation between ERa and CLOCK 
exists in these cells. Although the data appeared to indicate that 
ERa played a role in upregulating the expression of CLOCK in 
ERot-positive breast cancer cells, and that such regulation could be 
stimulated by E2, whether this mechanism is also important in 
normal breast cells needs to be addressed by further study. 

A possible model depicting how E2-ERot signaling is coupled to 
the machinery of circadian clock is shown in Figure 7. In this 
model, CLOCK, BMALl and PER2 transcription can be modulated 
by E2-ERa signaling. E2 enhances the sumoylation of CLOCK 
and the interaction of CLOCK with ERot. Sumoylated CLOCK 
may also function to increase the transcriptional activity of ERot. 
Meanwhile SENPl, which has been identified as a protease that 
desumoylates CLOCK may play a role in regulating the status of 
CLOCK sumoylation [27]. E2 also stimulates the transcription of 
PER2, leading to the accumulation of its transcript, and hence 
increasing the level of PER2 protein in the cytoplasm. PER2 will 



then be transported into the nucleus where it may inhibit the 
transcriptional activity of ERot [25]. Overexpression of PER2 has 
been shown to inhibit cell growth and the rise of clonogenic cells in 
breast cancer cells [25]. Moreover PER2 also inhibits the 
transactivation of the circadian proteins, CLOCK and BMALl. 
As the whole pathway is a network, the transcription of BAMLl is 
also regulated by E2 and ERa [63] . These previous studies seem to 
suggest that there is a closely relationship between circadian 
rhythm and E2-ERa signal pathway. The transcriptional regula- 
tion CLOCK mediated by E2-ERa signaling demonstrated in the 
current study would provide a positive contribution to the further 
understanding of the molecular mechanism by which E2 alters the 
circadian rhythm in behavior, physiology, and reproductive 
functions in mammals. At the same time, the crosstalk between 
E2-ERa signaling and CLOCK would add to the complexity of 
the mammalian circadian clock feedback loop. 

In conclusion, our results indicated that CLOCKk a downstream 
transcriptional target of ERa, and this provided potential insights 
into the connections between E2-ERa signaling and circadian 
rhythm, and showed that CLOCK may be an integral part of the 
series of genes that constitute the responsiveness of cells to the 
presence of estrogen, functioning as part of the network of 
transcriptional events governed by ERa. This will serve as a step 
forward in unraveling the complex mechanisms involved in the 
development of breast cancer involving a clock gene. 
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